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Despite observed strong correlations between central supermassive black holes 
(SMBHs) and star-formation in galactic nuclei, uncertainties exist in our un- 
derstanding of their coupling. We present observations of the ratio of heavily- 
obscured to unobscured quasars as a function of cosmic epoch up to z~3, and 
show that a simple physical model describing mergers of massive, gas-rich 
galaxies matches these observations. In the context of this model, every ob- 
scured and unobscured quasar represent two distinct phases that result from 
a massive galaxy merger event. Much of the mass growth of the SMBH occurs 
during the heavily-obscured phase. These observations provide additional ev- 
idence for a causal link between gas-rich galaxy mergers, accretion onto the 
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nuclear SMBH and coeval star formation. 

While unobscured quasars (i) have been known for a long time (2) and their statistical 
properties are well-studied (5), the numbers of the heavily obscured quasar population and its 
variation with cosmic epoch are still strongly debated. This population has been uncovered 
using multi-wavelength selection techniques that simultaneously exploit X-ray (4), optical (J) 
and mid-infrared (6) wavelengths. As a result of the efficiency of these techniques, the sample 
sizes of obscured quasars is growing substantially. The existence of a large number of heavily 
obscured quasars at 2;~2 was predicted by early AGN population synthesis models that success- 
fully explain the generation of the X-ray background (7). However, their space density cannot 
be constrained by these calculations (8). What population these obscured sources evolve into or 
proceed from is poorly understood at present. 

The link between Ultra-luminous Infrared Galaxies (ULIRGs) and quasars was first sug- 
gested by Sanders et al. (9). There is substantial observational evidence that ULIRGs, at least 
locally, are the product of the gas-rich merger of two massive (M^IO^^ Mq) galaxies {10). 
The merger process is believed to switch on accretion onto the central black hole as it provides 
efficient transport of gas to the nucleus {11). The gas funneled to the center is expected to fuel 
the supermassive black hole and induce star formation. The origin of the infrared luminosity 
of these sources, whether they are powered primarily by star formation processes {12) or the 
AGN (9) activity, is still debated. Here, we present recent measurements of the space density 
of heavily-obscured quasars as a function of redshift, and estimate the duration of the obscured 
stage by comparing models with observations. 

Mainly due to the effects of dust and gas obscuration at most wavelengths, finding heav- 
ily obscured quasars {13) is a challenging task which has prevented the identification of larger 
samples, in particular, at high redshifts. Measuring their space density is even more difficult, as 
it requires a good knowledge of the selection function and observational biases. We have com- 
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piled observations of obscured quasars selected at various wavelengths using different tech- 
niques, described in detail in the SOM, including spectral fitting in X-rays {14) and infrared 
selection (15-17). X-rays, specially at rest-frame energies greater than 10 keV, are not signif- 
icantly affected by obscuration. In addition, most of the absorbed energy is later re-emitted 
at infrared wavelengths. Thus, these techniques permit an estimate of the number of heavily- 
obscured quasars at 2;>1. In the local Universe, the space density of ULIRGs (18), which can 
be used as an indicator of the total number of quasars, implies that the ratio of heavily-obscured 
to unobscured quasars at z=0.l is ~1 (Fig. 1). 

Both in the local Universe and at 2;~1, heavily-obscured AGN are bright at mid- and far- 
infrared wavelengths and also potentially in the sub-millimeter but are optically faint (Fig. 2). 
At low redshift, an ULIRG is clearly the product of a galaxy merger, while at high redshift 
a merger is suggested, but deeper data are needed to confirm this hypothesis. The rest-frame 
optical images of six heavily-obscured quasar candidates show indications of ongoing major 
mergers and interactions (Fig. 3). 

In order to interpret these observations and the evolution of these populations, we started 
with the standard ansatz that the gas-rich major merger of two massive galaxies produces one 
newly-fueled quasar. This triggered quasar is originally obscured by the surrounding gas and 
dust (9), in some cases reaching Compton-thick levels, (i.e., where the optical depth for Comp- 
ton scattering is greater than one.) After a time At ~10^-10^ years (19), which we estimate 
independently below, most of the dust and gas are removed from the central region and the 
quasar becomes unobscured. 

To test this simple prescription, the calculated ratio of heavily obscured to unobscured 
sources from merger rates of massive gas-rich galaxies needs to match the observed ratio of 
obscured to unobscured quasars. This ratio can be calculated using the merger rate as a function 
of cosmic time in the context of the hierarchical cold dark matter (LCDM) structure forma- 
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tion paradigm. Using the assumptions described above, we estimated the ratio of obscured to 
unobscured sources as: 

AT ~ N f-y') ' 

-''unobsc ''unobsc I, ^ J 

where A^obsc is the space density of heavily-obscured quasars; d'^Mevgev/dNdt is the merger 
frequency per galaxy per unit time; A^gai and A^unobsc are the space densities of massive galaxies 
and unobscured quasars, respectively, and /g is the average fraction of gas-rich galaxies. The 
major merger frequency per galaxy in the LCDM paradigm can be parameterized as a power- 
law in (I-1-2;) with a mass-dependent exponent of ~ 1.5 (20). This form is derived from model 
parameters constrained by observations. In order to estimate the space density of galaxies above 
a threshold stellar mass, we used the median mass measured for ULIRGs found in the Cosmic 
Evolution (COSMOS) survey. The median mass increases with decreasing redshift, going from 
~10^^Mq at zr^2 to IO^^^Mq at z=O.S (Fig. SI). We then incorporated this limiting mass 
into the stellar mass function computed by Marchesini et al. (21) to obtain the space density 
of massive galaxies as a function of redshift. Rather than directly estimate the gas content 
of high redshift galaxies, which is currently observationally impossible at these redshifts, we 
used the average star formation rate as a proxy for the evolution of the fraction /g of gas-rich 
galaxies. The evolution of the star formation rate can be approximated as (1+2;)^ up to 2;~2, 
remaining mostly flat at higher redshifts, based on ultraviolet observations of galaxies up to 
2;~2.5 (22). Finally, the space density of unobscured quasars, A^unobsc(^) has been measured 
by both X-ray (23, 24) and optical (25) surveys, and consistent results are found with these two 
methods. These are all the ingredients needed to compute the expected fraction of obscured to 
unobscured quasars, wherein At in equation (1) can be determined as a free parameter. The 
redshift dependence of each of these components is shown in Fig. S2. 

The estimates from our simple scenario are consistent with the observations (Fig. 1), in par- 
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ticular considering the steep evolution in the relative number of obscured sources from 2;= 1.5 
to 3. This rapid increase is not predicted or expected from existing AGN luminosity functions 
(Fig. 1). That is, extrapolations from the behavior of less-obscured lower-luminosity sources 
(with observed column densities <10^^ cm~^) do not match current observations, in particular 
at z>2. This suggests the existence of a different channel for the triggering mechanism for 
quasars, compared to lower-luminosity AGN. The best-fit value of At we obtained is 96±23 
Myrs (90% confidence level). This is very similar to the current best estimates of quasar life- 
times in their optically bright, unobscured phases [10-100 Myrs (2(5)], indicating that these 
sources spend roughly half their life in the obscured phase. 

Having determined that quasars that are fueled by the merger of massive gas-rich galaxies 
spend comparable amounts of time in the the obscured and unobscured phases, we proceeded 
to estimate the implications for the mass accretion onto the nuclear supermassive black holes 
during these two stages. Assuming a typical accretion efficiency of £=0.1, the mass growth of 
the supermassive black hole due to accretion is given by: 



where T is the duration of the entire accretion episode that includes the obscured and un- 
obscured phases. The typical bolometric luminosities of these sources span the range 10^^- 
10^^ erg/s (-10^2.1014 l^)^ ^hile they have black hole masses of Mbh^IO^-IO^om© (27, 28). 
Given that the typical duration of the total (both obscured and unobscured) luminous quasar 
phase is T^lx 10^ years, it is possible for a quasar to build most or all of the black hole mass 
in a single event, which is triggered by a major merger as suggested here. 

In spite of copious accretion, due to their high redshifts and strong absorption, this unex- 
pected population of obscured quasars are not key contributors to the extragalactic X-ray back- 
ground radiation. Their total contribution is estimated to be ~l-2% {8). However, they com- 




(2) 
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prise a significant fraction of the super-massive black hole mass density in the Universe (29). 
Adding the extra obscured accretion reported here, which lasts as long as the optically bright 
phase, increases our original estimate of the integrated black hole mass density at z=0, by 4%, to 
4.5 X 10^ MgMpc"^ (8). Including this additional contribution, the integrated black hole growth 
in the obscured quasar phase is 1 .3 x lO^M^Mpc^'^, or ~30% of the total black hole mass den- 
sity at ^=0, in agreement with our conclusion that the obscured quasar phase can harbor a large 
fraction of the black hole growth (30). Our results are in agreement with recent estimates (26) 
that suggest an average accretion efficiency of <10% even accounting for heavily obscured 
accretion. 
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Redshift 

Fig. 1. The ratio of heavily-obscured to unobscured quasars as a function of redshift. 
Measurements of the space density of obscured quasars at high redshift were obtained from 

X-ray — green triangles (14) — and mid-IR imaging — blue pentagon (16) and black 
squares (17) — and spectroscopy — brown circle (15) — selection techniques. For the 2;^0 

measurement we used the luminosity function of local ULIRGs (18), assuming that each 
ULIRG is either a heavily-obscured or an unobscured quasar. The solid black line shows the 
heavily-obscured to unobscured quasar ratio expected from AGN luminosity functions derived 
from hard X-ray observations (32), while the red solid line corresponds to the ratio obtained if 
every gas-rich major merger of two massive galaxies generates a heavily-obscured quasar, 
which after a time At^96 Myrs becomes unobscured. Dashed lines show the uncertainty in 

this relation, at the 90% confidence level. 
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Fig. 2. Optical images of two examples of heavily-obscured, luminous AGN observed with the 
Hubble Space Telescope. Filters used are F435W (B) and F814W (I) and V,I,z in the left and 
right panels respectively, while the physical size of both images is ~90 kpcx90 kpc. Contours 

in the right panel show the Spitzer emission. The left panel shows the ^=0.03 prototype 
ULIRG Mrk 273, while in the right panel we present a high redshift heavily-obscured quasar 
candidate in the GOODS-South field, with no X-ray detection. The former is part of the IRAS 
1 Jy sample (18) and has clear indications of a major merger, while X-ray observations with 
Suzaku reveal Compton-thick obscuration levels (A^h>10^'^ cm~^) and intrinsic quasar-like 
luminosities (33). The high-z source was selected as a Compton-thick AGN candidate based 

on its high mid-IR to optical flux ratio and red optical/near- IR colors. Left panel credit: 
NASA, ESA, the Hubble Heritage Team (STScI/AURA)-ESA/Hubble CoUaboration and A. 
Evans (University of Virginia, Charlottesville/NRAO/Stony Brook University). 
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Fig. 3 Rest-frame optical images of 6 mid-lR- selected heavily-obscured quasars at in the 
ECDF-S region. Top panel show images obtained using the HSTAVFC3 camera using the F, J 

and H observations of the Ultra-Deep (left image) and GOODS fields. The bottom images 
were made combining data in the i?, J, and K bands obtained from ground-based telescopes, 

hence with ~10x worse spatial resolution than the HST images. All images are 15"x 15". 
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Observations of obscured quasars 

Finding obscured sources is a difficult task. So difficult that in fact for a while it was suspected 
that these sources did not exist. However, now mostly due to extensive multiwavelength surveys 
a significant number of these sources have been identified. In this work, we have compiled 
observations of obscured quasars selected at various wavelengths using different techniques, 
which are described below. 

In the nearby Universe, we can estimate the number of heavily-obscured quasars from the 
space density of ULIRGs. In particular, we base our calculation in the IRAS Revised Bright 
Galaxy Sample (S15). This sample includes both obscured and unobscured quasars, since only 
an IR selection was done. A random subsample of 16 sources with Lir>10^'^Lq was observed 
in X-rays by Chandra (S9). Of those, 9 have clear AGN signatures, namely a flat X-ray spec- 
trum, and hence the fraction of AGN in this ULIRG sample is at least 56%. However, for the 
remaining 7 sources a strong Fe K emission lines is detected in the stacked X-ray signal. Since 
detection of Fe K in emission is considered an unambiguous AGN signature, we can conclude 
that the fraction of sources containing an AGN, even in the sample of 7 sources not dominated 
by the AGN emission, should be very high. Similarly, although at higher redshifts, from a sam- 
ple of 70 //m-selected sources in the Cosmic Evolution Survey (COSMOS) field detected with 
Spitzer, J. Kartaltepe et al. (Sll) conclude that the fraction of sources with AGN signatures at 
LiR> IO^^Lq is higher than 80% and 100% for Lir> IO^^Lq. Hence, our assumption that most 
(if not all) of the ULIRGs contain an AGN is well justified. 

Heavily-obscured quasars can in principle be selected by their X-ray emission. The effects 
of obscuration are more important at softer energies, while harder X-ray photons can in principle 
escape more easily. This is why it is possible to study even Compton-thick sources at high red- 
shifts using the Chandra bands, which trace rest-frame energies > 10 keV at z> 1 . There are two 
types of Compton-thick quasars based on their X-ray spectra: transmission-dominated sources, 
characterized by a photoelectric cutoff at £'~10 keV corresponding to a cm~^ and 

in which the intrinsic emission is visible at higher energies, and reflection-dominated quasars, 
which are in general more obscured sources in which the intrinsic emission at X-ray wave- 
lengths is completely absorbed and only a weak (~l-2% of the intrinsic luminosity) reflection 
component can be detected. 

Using the deep Chandra observations available in the Chandra Deep Field South (CDF-S), 
R Tozzi et al. (S17) identified 10 Compton-thick candidates (both transmission and reflection 
dominated) with intrinsic luminosity Lx>10^'^ erg/s in the 0.73<z<4.29 range. These sources 
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were used to compute the space density of Compton-thick AGN in two redshift bins at and 
zr^l, as shown in Fig. 1 in the main text. 

As most of the absorbed energy is re-radiated at infrared wavelength, obscured AGN are in 
general very bright in the mid-IR bands (e.g., SI 8). Using a combination of optical and mid-IR 
spectroscopy, together with X-ray data, seven Compton-thick quasar candidates at z~2 were 
identified {SI). The space density of Compton-thick quasars was estimated by these authors 
from the four sources found in the GOODS-North field to be ~0.7-2.5 x 10~^ Mpc~^. 

Using a method based only on optical to mid-IR photometry, F. Fiore et al. {S3) showed 
that the sources with a high mid-IR to optical flux ratio present at the same time significant X- 
ray emission, higher than the expected value for pure star-forming galaxies and consistent with 
Compton-thick obscuration levels. This selection method was then applied to the sources in the 
COSMOS field {S4). This provides a space density of Compton-thick quasars of 4.8±1.1 x 10~^ 
Mpc~^ in the 1.2<z<2.2 redshift bin. As can be seen in Fig. 1 of the main paper, this value is 
in very good agreement with the density derived from X-ray observations {S17) for sources at 
similar redshifts using a completely independent selection technique. 

The same optical to mid-IR selection technique was used in the Extended CDF-S (ECDF-S) 
field {S19), which has very deep Spitzer and Chandra imaging over a ~0.3 deg^ region of the 
sky. A total of 211 Compton-thick quasar candidates were found using the selection scheme 
of F. Fiore et al. {S3). Thanks to the larger number of sources in this sample, it is possible to 
measure the density of heavily-obscured sources in four redshift bins, from 2;~1.5 to Z'^1.9. 

Estimating the ratio of obscured to unobscured quasars from 
galaxy mergers 

In order to calculate the redshift dependence of the ratio of obscured to unobscured quasars 
we start with the following reasonable and justifiable assumptions. New quasars are ignited by 
the gas-rich merger of two massive galaxies. The product of this merger produces a ULIRG, 
which harbors at its center a heavily-obscured quasar. After a time At most of the surrounding 
dust and gas are blown out by the radiation pressure of the rapid on-going accretion, after 
which an unobscured optically bright quasar is visible. In this scenario, the ratio of obscured 
to unobscured quasars as a function of redshift is simply given by the merger rate of gas-rich 
galaxies (equation 1) in the main text of the paper. Here, we describe and evaluate each term in 
that equation and its redshift dependence in detail. 

One of the fundamental ingredients in this calculation is the evolution of the merger rate 
and its dependence on other parameters like e.g., galaxy mass. This is still a controversial topic. 
While work based on galaxy morphology to identify mergers reported no redshift dependence 
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(e.g., S12), others based on galaxy pairs (e.g., SIO) found a strong redshift dependence of the 
form ~(l+z)^. A compilation by P. Hopkins et al. (S7), which includes the sources in the 
two samples mentioned before, concluded that once the differences in median galaxy mass and 
mass ratio (minor versus major mergers) are accounted for, the dependence of the merger rate 
on redshift is confirmed at a high statistical significance. Following this work, the major merger 
rate per galaxy as derived from simulations constrained by observations is parametrized as: 



where 



- MMuMi + zf'""-^ (3) 
X(MM,J = 0.02[l+(^^3^nGyr-' 

and 

M 

/3(Mm.) = 1.65- 0.15 log(^-f5j5^). 

Uncertainties in these parameters are a factor of ~2 in A(MMin) and 0.2 in ^(MMm)- The 
minimum stellar mass in our definition of "massive" galaxy is a function of redshift, as shown 
in Fig. SI, and can be parametrized as 

MMin(^) = 5 X 10^1(1 + z)-'-'Mq. (4) 

Hence, A{M) ranges from 0.12 Gyrs~^ at z=0 to 0.06 Gyrs~^ at z=3, while /3{M) goes from 
1.44 at z=0 to 1.58 at z=3. The resulting redshift dependence of the merger rate per galaxy is 
shown in Fig. S2. 

The space density of galaxies as a function of stellar mass and redshift was determined 
observationally recently by D. Marchesini et al. (S13). They report that the stellar mass function 
is well-fit by a Schechter (S16) function up to 2;~4 in which the normalization is a function 
of redshift. To calculate the merging rates of massive galaxies, we use the limiting stellar 
mass given by equation (2). The resulting space density of massive galaxies as a function of 
redshift is shown by the blue line on Fig. S2. In order to estimate the number of gas-rich 
galaxies the extinction-corrected star formation rate density as a function of redshift is used 
as a proxy of available gas mass at each epoch. This is done because direct measurements 
of gas contents, especially in high redshift galaxies, are very difficult and hence we have to 
rely on secondary indicators. In particular, we use the measurements obtained from the Great 
Observatories Origins Deep Survey (GOODS) observations (52). A good fit to these data, 
normalized to the high-redshift value, is obtained using the following expression 

Gas rich(z) = S ^ z>2 ' 
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This dependence (shown in Fig. S2) is consistent with results obtained from hydrodynam- 
ical simulations (S6). Finally, the number of unobscured quasars is obtained directly from 
observations, in particular from X-ray luminosity functions (S5). Since this selection was made 
at soft X-ray energies, this ensures that only unobscured sources are included. Only sources 
with Lx>^0'^'^ erg/s, corresponding to Lboi>10''^ erg/s, were considered. If instead an optical 
luminosity function (e.g., S14) is used, assuming in this case an equivalent magnitude cut of 
Mi<-23 corresponding to a similar luminosity threshold, the results are consistent. As can be 
seen in Fig. S2, the quasar number density has a strong peak at 2;~1.5, followed by a shallow 
decline towards higher redshifts. 

The redshift dependence of each of the components in our calculation is shown in Fig. S2. 
Both the number of unobscured quasars and massive galaxies show steep declines from 2;~1 to 
z=0. Meanwhile, the merger rate increases steadily by a factor of ~3 from z=0 to z=3. 
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Fig. SI. Median stellar mass as a function of redshift for the host galaxies of the ULlRGs 
(LiR>lO^^L0) in the COSMOS field. Stellar masses were derived by performing spectral 
fitting of the observed-frame optical and near-IR photometric data (S8). Vertical error bars 
were obtained from the dispersion in each bin, while the horizontal error bar show the width of 
each bin. The dashed line shows the best-fitting power-law, 5 x 10^^(1-1-2;)^^^ Mq, which is 
used to compute the number of galaxies hosting ULIRGs as a function of redshift. 
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Fig. S2. Evolution of the different ingredients used in our calculation of the expected ratio of 
obscured to unobscured quasars as a function of redshift. All curves were normalized to their 

z=0 value. Black line shows the evolution of the density of unobscured quasars, while the 
redshift dependence of the major merger frequency per galaxy, density of massive galaxies and 
fraction of gas-rich galaxies is shown by the red, blue and green lines. 
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